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ABSTRACT 

We continue our earlier studies of quasi- periodic oscillations (QPOs) in the 
power spectra of accreting, rapidly-rotating black holes that originate from the 
geometric “light echoes" of X-ray flares occurring within the black hole ergo- 
sphere. Our present work extends our previous treatment to three-dimensional 
photon emission and orbits to allow for arbitrary latitudes in the positions of the 
distant observers and the X-ray sources in place of the mainly equatorial posi- 
tions and photon orbits of the earlier consideration. Following the trajectories of 
a large number of photons we calculate the response functions of a given geome- 
try and use them to produce model light curves which we subsequently analyze 
to compute their power spectra and autocorrelation functions. In the case of an 
optically-thin environment, relevant to adveetion-dominated accretion flows, we 
consistently find QPOs at frequencies of order of ~ kHz for stellar-mass black 
hole candidates while order of ~ mHz for typical active galactic nuclei (~ 10 7 i\/r,) 
for a wide range of viewing angles (30° to 80°) from X-ray sources predominantly 
concentrated toward the equator within the ergosphere. As in our previous treat- 
ment, here too, the QPO signal is produced by the frame-dragging of the photons 
by the rapidly-rotating black hole, which results in photon “bunches'’ separated 
bv constant time-lags, the result of multiple photon orbits around the hole. Our 
model predicts for various source/observer configurations the robust presence of 
a new class of QPOs, which is inevitably generic to curved spacetime structure 
in rotating black hole systems. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: 
galaxies — stars: oscillations 
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1. Introduction 

One of the more interesting findings of the X-ray timing analyses of the light curves of 
accreting compact objects in the past two decades has been the discovery of quasi-periodic 
oscillations (QPOs). These are broad features in the power spectra of these sources that 
range up to kHz frequencies in neutron star low-mass X-ray binaries (see, e.g., van der 
Klis 2000) and up to tens and even hundreds of Hz in accreting black hole binary systems 
(see, e.g., Cui et al. 1998 and Strohmayer 2001a. b), including QPO pairs with 2:3 frequency 
commensurability in the power spectra of certain black hole candidates (e.g. XTE J1550-564 
and GRO J1655-40). Besides galactic binary sources, the presence of QPOs has also been 
reported in at least one ultra-luminous X-ray sources (ULX), namely NGC 5408 X-l (see 
Strohmayer et al. 2007) at v ~ 20 mHz; in the context of active galactic nuclei (AGNs) 
Gierlinski et al. (2008) recently discovered a ~ 1 hour X-ray periodicity in the narrow-line 
Seyfert 1 (NLS1) RE J 1034+396, in which the periodicity is apparent in its light curve in 
distinction with those of the galactic binary black hole systems mentioned above that do not 
exhibit any such behavior. 

The nature of QPO phenomenon in accreting black hole systems is at present largely 
unknown; it is also possible that this phenomenon is not the manifestation of a unique 
process, but that the observed diversity of QPO frequencies and sites may be due to a 
multitude of processes specific to the particular frequency and site. For the case of QPOs 
associated with black hole systems, the absence of an underlying, rotating solid surface 
object, implies that plausible explanations involve by necessity processes associated with 
the surrounding accretion disks. These include, among others, accretion disk precession due 
to the frame-dragging effects of a rapidly-rotating black hole (e.g. Cui et al. 1998; Merloni 
et al. 1999; Schnittman et al. 2006) (also see, e.g., Stella & Vietri 1998, for the same effect 
invoked to account for the kHz QPOs of Low Mass X-ray Binaries), accretion disk oscillatory 
modes (e.g. Nowak et al. 1997; Abramowicz & Kluzniak 2001; Kato 2001; Donmez 2007, for 
diskoseismology), or the distribution of X-ray emitting “blobs” spanning a limited range 
around some specific radius of the Keplerian disk surrounding the black hole (e.g. Karas 
1999). There have also been studies of the correlations of the QPO frequencies with the 
properties of the associated X-ray spectra (e.g. Titarchuk fc Fiorito 2004). 

An altogether different notion to QPO origin has been put forward recently by Fuku- 
mura A Kazanas (2008, hereafter FK08). These authors proposed that high frequency QPO 
(HFQPO) could be produced as a consequence of light, echoes of X-ray flares occurring within 
the ergosphere of a Kerr black hole: they showed that due to frame-dragging of individual 
flare photons, a finite number of them reach far away observers after an integer number of 
additional orbits around the black hole to produce a geometry-induced “light echo” of the 
original flare (see. also. Meyer et al. 2006: Bursa et al. 2007, for a similar discussion). Study- 
ing two-dimensional (2D) near-equatorial photon propagation in a fully general relativistic 
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calculation, FK08 showed that dragging of inertial frames leads approximately 15% of the 
photons to a specific direction at infinity (an observer) after a time lag of At ~ 14 M (where 
M is black hole mass), independent of the relative position between the observer and the 
photon source. The independence of the lag on the relative phase between the observer and 
the source then guarantees a second peak in the Autocorrelation Function (ACF) at lag 
r ~ 14 M, which according to Fourier analysis, results in a prominent HFQPO signal at a 
frequency t'qpo — 1/ At ~ 1.4 (10M©/M) kHz (where M© is the solar mass), even for a light 
curve that consists of flares randomly distributed within the ergosphere. Since the closest 
an accretion disk can approach the black hole is the radius of Innermost Stable Circular 
Orbit (ISCO), this effect is possible only when the ISCO lies within the black hole ergo- 
sphere, a condition that limits this effect to black holes with dimensionless spin parameter 
a/M > 0.94. 

The study of FK08 showed the prominent role of frame-dragging in producing lags 
independent of the relative source - observer position; as shown there, the absence of frame 
dragging for flares in a Schwarzschild geometry leads to lags which depend on the relative 
observer - source position and, for random flare positions, to absence of QPO features in the 
corresponding power spectra. 

The original study of FK08 outlined the fundamental aspects behind this work, namely 
the importance of frame-dragging and the persistence of QPO presence, even for random 
positions of the X-ray flares within the ergosphere. However, it was restricted in its scope 
in that it was constrained to observers and sources near or on the black hole equatorial 
plane, i.e. 6 S ~ 9 a ~ 7t/2. Realistic constraints demand an enlarged study to include, at 
a minimum, observers at larger latitude positions, not least because it is expected that the 
column density of sources for lines of sight near the equator are sufficiently high to preclude 
propagation of photons without interaction with the surrounding matter; this could then 
lead either to their absorption or to their scattering and introduction of additional, random 
lags to their trajectories erasing the otherwise potential QPO signal. The goal of the present 
work is to remedy this deficiency by expanding the study of the corresponding lags to three- 
dimensional (3D) geometries allowing thus for arbitrary latitudes of both the observers as 
well as of the source of photons. 

We would like to stress that the present models of “geometry-induced” QPOs (and also 
those of FK08) are not meant as an account of the so-called HFQPOs presented in the 
literature to date (e.g. Cui et al. 1998; Strohmaver 2001a,b). For one thing, the predicted 
QPO frequencies are 5—10 times higher than the observed frequencies mentioned above 
(for the same black hole mass), and. for another, they do not produce their 2:3 frequency 
commensurabilitv. Instead, the predicted QPOs should be viewed as an altogether new 
feature indicating the intrusion of an accreting gas within a black hole ergosphere, thereby 
providing an unequivocal evidence for the presence of rapidly rotating black holes. 
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Our paper is organized as follows: In §2 we provide a general description of our 3D 
QPO model, the details of the photon kinematics and response function of the system, 
and a prescription for constructing stochastic model light curves. In §3 we compute the 
autocorrelation functions (ACFs) and power spectral densities (PSDs) of the model light 
curves to show that they generally exhibit the QPO features as anticipated. Finally, in §4 
we review our results, make contact with observations along with the limitations of our model, 
and conclude by discussing prospects of future work. Finally, in the Appendices we discuss 
the details the notion of locally isotropic emission in Kerr geometry and provide explicit 
functional forms relating the photon impact parameters to their local emission angles. 


2. Model Details 

The present treatment follows closely that of FK08, but generalizing the photon tra- 
jectories to 3D, thus allowing for positions of both the source and the observer at latitudes 
other than near equatorial. So, we consider photon emission in Kerr geometry of dimension- 
less angular momentum of a black hole a, described in Boyer- Lindquist coordinates (r, 0,4>). 
As in FK08, here too we employ the usual geometrized units (G = c = 1 where G is the 
gravitational constant and c is the speed of light) and normalize distance and time with the 
black hole mass M. For example, 1 M corresponds to 1.5 x lO 5 (M/A/ 0 ) cm in distance and 
5 x 10 ~ 6 (M/M 0 ) sec in time. We consider the observer position at an arbitrary location 
( r 0 ,9 0 ,(f ) 0 ), while correspondingly denoting the source coordinates as (r s ,9 s ,(/) s ). This ar- 
rangement, therefore, allows for source locations not only at random azimuths but also at 
random polar angles. In thin accretion disks, the flares are likely confined to small latitudes, 
i.e. 6 S ~ 7 t /2, but for quasi-spherical accreting flows such as advection-dominated accretion 
flows (ADAFs; e.g. Narayan & Yi 1994), the overall source latitude can be much higher. As 
discussed in following subsections, in this last case, one can, in addition, weigh the emission 
at each latitude to fit one's preferred flow model. 


2.1. X-ray Flares 

The X-ray emission of the sources we consider is regarded as the incoherent sum of 
a large number of instantaneous (i.e. short compared to the dynamical time) flares. We 
assume the typical duration of such a flare to be ~ 1M. presumably originating from local 
magnetic reconnections (e.g. Svensson & Zdziarski 1994) or standing shocks in relativistic 
magnetized accretion (e.g. Takahashi et al. 2002: Fukumura et al. 2007). Concerning gen- 
eration of X-ray flares within the ergosphere, Karas & Kopacek (2008) have attributed it 
to magnetic reconnection induced by frame-dragging, while Koide & Arai (2008) have con- 
sidered magnetic reconnection within the ergosphere as a means of extracting energy from 
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Fig. 1. — (a) Schematic view of the system considered in this paper, (b) Poloidal projection 
of the geometrical setup in this problem. 

the black hole rotation. In general, these X-ray flares do not necessarily have to be confined 
exactly to the equatorial plane, and for this reason, in this work, we relax this restriction 
on the source vertical location. Therefore, in the specific cases that the emitting plasma is 
optically thin, the observed signal wall contain contributions from sources in both upper and 
lower hemispheres. 

It is important to mention that in this model both polar and azimuthal positions (9 a , 6 S ) 
of each X-ray flash are randomly assumed (as long as they lie within or close to the ergo- 
sphere) while their cylindrical radius r s sin 6 S is fixed for simplicity to the radius of the ISCO 
of a Kerr hole wdth a/M ~ 0.99, i.e. to r s sin 0, ~ 1.455M. One of our primary interests 
is thus to probe the dependence on the source and the observer latitudes as illustrated in 
Figure lb. To avoid the introduction of artificial coherent signals associated with orbital 
motion of these emitting sources, we also consider the flares to be randomly distributed in 
time in the following fashion: each X-ray burst occurs with a Poisson distribution of mean 
value equal to the Keplerian orbital timescale; AT = T orb | ln{rnd(0, 1)}| where T mh is the 
orbital time of the source at ( r s ,6 s ) and rnd(O.l) is a random number between 0 and 1. We 
define that Torb = 27r{(r s sin 0 S ) 3 / 2 + aM 1 ' 2 }/M 1/2 . 

The orbiting X-ray emitting plasma is assumed to be in Keplerian motion, either below 
or above an equatorial plane at 0 s i where i denotes the i-th source, i.e. rotating with the 
local Keplerian angular velocity (modified to take into account the black hole spin a) 

. I / 1 /' 2 

(rsiii0) 3/2 + aM 1 / 2 

with all the quantities being evaluated at the source position. In the lab frame or the locally 
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non-rotating reference frame (LNRF) each source has a total three- velocity of 


v a 


A sin 6 
EA 1 / 2 


(fi — id) 


J source 


j 


( 2 ) 


where we define A = r 2 — 2 Mr + a 2 , E = r 2 + a 2 cos 2 0, A = (r 2 + a 2 ) 2 — a 2 A sin 2 9, and 
frame-dragging of the local inertial frame is denoted by u = 2Mra/A. In the LNRF, as one 
can see, rotation of a local inertial frame lu has been subtracted by definition. 


2.2. Null Geodesics 


Neglecting external interactions between photons and particles, each null geodesic (i.e. 
photon orbit) is uniquely characterized by two constants of motion, £ and r/, where £ is 
the axial component of angular momentum and r/ closely related to its polar component 
(e.g. Bardeen et al. 1972; Chandrasekhar 1983). For given r/ and £, a photon orbit in Kerr 
geometry is generally governed by the following equations of motion (e.g. Bardeen et al. 
1972; Chandrasekhar 1983); 


(r 2 + a 2 )E + aA£ — a 2 A sin 2 # 
= _ 

r 2 E 2 = H 2 -A{r, + (£-a) 2 }, 

#E 2 = ?/ + (a 2 — £ 2 esc 2 9) cos 2 9 , 

■ a(2Mr — a£) + £A esc 2 9 


(3) 

(4) 

(5) 

( 6 ) 


where E = r 2 +a 2 — a£ and an overdot denotes derivative with respect to the affine parameter. 
Note that equations (4) and (5) provide only the squares of r and 9. In order to avoid the 
issue of determination of appropriate signs for r and 9, for the actual numerical determination 
of the orbits we use instead the second derivatives of equations (4) and (5) 




9 = 


2 rrr, — r (r — M) {>/ + (£ — a) 2 } — E |2rr — a 2 sin(2#)/9 j r 2 

— 

£ 2 cot 9 esc 2 99 — a 2 sin(29)9/2 — E |2rr — a 2 sin(2fl)ti| 9 2 


(7) 

( 8 ) 


The initial conditions necessary for the integration of the orbit equations involve, besides 
the source position. (r s .9 s ), also the two angles of photon directions and the values of r and 
9: the last four quantities are not independent but they can be expressed in terms of (r,.. 9 S ) 
and the values of the impact parameters ij and £ as discussed in §2.3. 
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Fig. 2. — Geometry of local emission in the LNRF described in the text. 


2.3. Locally Isotropic Emission 

In the absence of any information concerning the angular distribution of the photons 
emitted at a given flare, the most conservative assumption is that the emission is isotropic 
in the rotating fluid rest frame. However, the proximity of the X-ray sources to the black 
hole and the associated source motion (u s ~ 0.6 c) and strong field geometry along with the 
fact that the coordinate system for the integration of the photon orbit equations is that of 
the LNRF, require that the notion of photon isotropy be transformed to this frame and that 
the impact of the Kerr metric on the notion of isotropy be considered with some care. Below 
we discuss our prescription of local isotropic emission, while more details can be found in 
Appendix A. 


The geometry of photon emission is depicted in Figure 2. We assume the black hole 
spin axis to be in the z-direction while the local fluid velocity v s in the y-direction; ip is the 
angle between the emitted photon and the instantaneous direction of source motion (y-axis) 
in the fluid frame, while \p' is the corresponding angle in the LNRF. For a given 'ip the precise 
direction of emission is determined by the angle \ confined in the plane orthogonal to the 
y-axis. The relation between -ip and x p' is then (the angle \ remains invariant) 


v) 


cos 


-1 


3 + cos ip 
1 i- 3 cos r 


( 9 ) 


where 3 = r s /<\ while the differential solid angle and opening angle are transformed respec- 
tively as 

1 — i' (1 — 

r/(cosr ) = — 3 —(/(cost’) . chi' ( 3 . v) ~ - dv . (10) 

v (l + jcosc') 2 v 1+3 cos ip v ' 


Isotropy in the fluid frame implies equal number of photons per unit solid angle .\ 0 . i.e. 
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Fig. 3. — Comparison between the fluid rest-frame and the LNRF with /3 = 0.5 for (a) 
Differential photon number distribution per emitting angle as a function of respective emit- 
ting angle in the fluid frame (dark) and the LNRF (gray), (b) Corresponding bin-corrected 
photon distribution (i.e. weighting factor) in the fluid frame (dark) and the LNRF (gray). 
We take Nq = 1000 and A ip = 7t/50. 


dN/dfl = Nq over the entire sky in the fluid rest frame. However, when computing photon 
geodesics we choose incrementally not the solid angle but the angles i p and y. Then, the 
number of photons emitted per increment of the polar angle ip is dN / dip — 2ttN 0 sin ip and 
therefore the increment of the angle x should be equal to 1/ (2txNq sin ip) so that each solid 
angle element at a given ip, x contains one photon. 

The Lorentz boost to the LNRF skews the photon distribution with angle resulting 
in larger number of photons being emitted along the direction of source motion with the 
distribution as a function of the polar angle ip (or ip') given by the function 


/(Art 


(1 — 3 2 ) sin ip' 
(1 — 3 cos ip') 2 


(1 + 3 cos C) sin^/ 

(i-V ) i/2 


mi o ■ 


(ii) 


In Figure 3a the photon distribution as a function of the polar angle ip for 3 = 0 (dark dots 
for the fluid frame) and also as a function of ■ ip ' for 3 = 0.5 (gray dots for the LNRF) is 
shown. We take .Y„ = 1000 and A tp = tt/o0. However, as indicated above in equation (10), 
the Lorentz boost transforms also the width of the the angular bin Ar: because of that, 
as shown in Appendix A. the number of photons per angular bin is invariant (as it should 
for photon number conservation), proportional to sin y AC. So each angular 1/in Ac' in 
the LNRF has the same number of photons as the corresponding bin An in the fluid rest 
frame, but these photons are now received at a smaller angle yP, given in terms of i< by 
equation (9). This distribution is shown in Figure 3b. where the number of photons per 
angular bin f(3 , p')Ay' = (dN /dv')Ay’' — (dN/d4')^y: •- for the same values of A' 0 and A< • 
as in Figure 3a - is shown as a function of the angle -tp (dark dots for the fluid frame) and 
(A (gray dots for the LNRF). It is apparent that a given number of photons per bin is found 




at smaller angles in the LNRF rather than in the fluid rest frame. 

For our photon orbit calculations we transform an equally spaced (discrete) set of angles 
ipi in the source rest frame to the corresponding angles ip'i in the LNRF and release at each 
such bin a number of photons equal to 2ttNq sin ip t . These photons are then distributed each 
in a single of as many equal intervals in the \ angle. The next step involves expressing local 
emitting angles S' and 7' (both measured in the LNRF) used in the integration of the photon 
trajectories in terms of the angles ij’' and y. The relations between these sets of angles can 
be obtained from trigonometric relations of the geometry shown in Figure 2 as 

cos ip' = sin S' sin 7' , cos<5 ; = sin ip' cos 7 , (12) 

where the first of these equations was obtained from a spherical triangle for the cosine of 
the angle ip' while the second one from relations of the projection of the geometry of the 
same figure onto the plane perpendicular to the y-axis (i.e. in the shaded ellipse of the same 
figure). 

The next step in the calculation takes into account the effects of the local geometry 
curvature in the definition of local isotropy (see Fukumura & Kazanas 2007, as an example 
of this issue in the simpler case of the Schwarzschild geometry) by relating the angles S' 
and 7' to the photon four-velocity components. In the LNRF the photon’s local emitting 
polar and azimuthal angles (S', 7') are related to the four-velocity components through the 
relations 

r 

v r = — — - = sin S' cos 7' , (13) 


AY' Z 6 


A t 


A sin# ( <b 


EA !/2 l t 


uj = sin S' sin 7' 


which lead to 


A \ 1/2 6 
A t 


— - AVW*smeo-<J ’ V ; 

evaluated at the emitting location (r s ,# s ,0). Using equations (3) through (6) we can elim- 
inate t. r. 0 and o from the above equations (16) and (17) in favor of the two impact 
parameters // arid £. We are thus led to two relations that connect the photon emission 
angles (S'. 7 ') to the photon impact parameters (£. //). In fact, one can invert these relations 
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Fig. 4. — (a) Locally isotropic photon trajectories from a source at (sin 9 S , (f> s ) = (89°, 0°) for 
an opening angle of 0 = 80° (ip' — 46°). We set a/M = 0.99 and N 0 = 100. (b) Observed 
(cumulative) photon distribution Af to t(^)/A 7 tot(90°) as a function of covering angle measured 
from the pole for various source geometries. It is normalized by the total photon number 
over 7t/2. We assume an equatorial source of 6 S = 89° (dark curve) and mid-latitude source 
of 53° (gray curve), in comparison with stationary flat-space case (dotted curve). 

to obtain (£,77) as a function of these angles, i.e. £ = £(<5',7 7 ) and rj = rj(5', 7') as shown in 
Appendix B. This now completes the problem of the choice of the orbit parameters consis- 
tent with isotropic emission in the fluid frame: For a given set of equally spaced values of ipi 
(locally isotropic source in the fluid rest frame) one can compute ip( and also determine the 
number of photons to be emitted and values of the azimuthal angle bins \, ; these are then 
used to compute the angles S' and 7', and then the corresponding impact parameters rj and 
£ along with the initial values of r(5', 7') and 6(8'. 7') that are used in the integration of the 
photon orbit by equations (7) and (8). 

An example of 3D photon trajectories from a locally isotropic emission is shown in 
Figure 4a where we consider emission with a half-opening angle of v = 80° (ip 1 = 46°) 
from a source at (6 s ,q s ) = (89°, 0°). Throughout this paper we assume a/M = 0.99 as in 
Figure 4a. unless otherwise stated. Recall that our model considers now X-ray sources from 
both the upper and lower hemispheres with respect to the equator because of the opticallv- 
thin assumption of surrounding matter (see Fig. 1). Hence, the net observed emission will 
consist of a sum of these two components. We have confirmed that these photon trajectories 
are indeed almost symmetric with respect to the equator as expected. It should be noted, 
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however, that for some particular emission angles this symmetry in trajectory with respect to 
the equator can be largely broken even for a small deviation of the source position from the 
equator. Due to the high velocity of the source (J3 ~ 0.6) for the values of r s and a chosen, 
most of the emission is beamed forward in the direction of instantaneous source velocity. In 
fact we find that many backward-emitted photons (i.e. of 0 > 90°) subject to the Doppler 
and frame-dragging effects of the fast black hole rotation are forced to corotate in the hole 
direction (see FK08 for details in 2D cases). We also find that backward-emitted photons 
with 0 > 140° end up crossing the the horizon, making no contribution to the observed 
signal. We should mention that the photon trajectories considered here also include many 
multiple orbits and higher-order photons; i.e. we also collect those photons that orbit around 
a black hole multiple times with some time-lags. The physical significance of these time-lags 
will be further discussed in §2.4. 

To illustrate the effects of Doppler beaming and the black hole’s strong gravity for 
the considered source location (r s sin 9 S — 1.45 A/), we present in Figure 4b the cumulative, 
normalized distribution of photons A r tot (D) /TV tot (9f) c ) received by a distant observer at r/M = 
100, as a function of the observer’s inclination angle 9„. In the figure the solid dark curve 
corresponds to a nearly equatorial source (6 S = 89°), the gray curve to a source at the edge 
of the ergosphere (0 S ~ 53°) for the fixed cylindrical radius of the source, while the dotted 
curve to the flat-space distribution. It is apparent in this figure that photons emitted from 
equatorial position (in particular) are directed primarily at low latitude directions. While in 
flat-space 50% of the photons (horizontal dotted line) are emitted between angles 0° — 60° 
as expected, in the equatorial source case of 9 S = 89° the half-way point has been extended 
to 9 ~ 82°. Similarly, while in flat-space polar angles 0° — 45° include ~ 30% of the entire 
photons emitted, in the present Kerr geometry they include roughly 10 times less. This is 
an important issue related to the equivalent width of the (broad) fluorescence Fe lines whose 
red wing is presumably emitted by plasma closest to the ISCO of fast-rotating black holes 
with spin parameter similar to that used in producing our Figure 4b. 


2.4. Response Functions and Light Curves 

Following a methodology similar to that discussed in FK08, we first compute a source’s 
response function seen by a distant observer whose radial and azimuthal positions are set 
in this work to (r 0 /M,o 0 ) = (100,0). This is a reasonable choice because at this distance 
the photon orbits are already nearly straight and thus their final (polar and azimuthal) 
directions are almost the same as those in flat spacetime. For a given source position at 
{O.s.h Os.i) around a rapidly-rotating black hole we collect isotropically emitted photons (in 
the local rest-frame of an orbiting source) keeping track of its arrival time t and final position 
(9. o) for each photon. As a result we can construct a response function seen by the observer 
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Time (in M) Time (in M) 

Fig. 5. — Synthetic response functions from a single X-ray source at 9 S = 89° with (a) 
(f> s ~ 4° and (b) ~ 180° seen by an observer of 0 o = 60°. Arrows indicate a constant lag of 
~ 144/. Values of the other parameters are the same as in Figure 4. 

Ij(t, 9 0 : 9 s ,i, which tells us photon counts as a function of photon’s arrival time from the 
i-tli X-ray source at (6 a ^, <p SJ ). For an optically-thin environment we have contributions from 
sources at different source latitudes for the same source azimuth, however, for simplicity we 
assume that the sources closer to the equatorial plane are dominant, although it is found 
that this assumption does not make much difference to our final results shown §3 below. 

The temporal resolution in the response function is taken to be \M (corresponding to 
5 x 10 5 sec for a 1OM 0 black hole), therefore the signal received by the observer is uniformly 
binned with this numerical accuracy. For computational purposes we approximate this signal, 
as precisely as possible, by a set of narrow Gaussians, also of width 14/. Note the response 
function from each single source contains a characteristic temporal profile due to distinct 
photon trajectories between the source and the observer, and the profile also strongly depends 
on the observer’s polar position 9 0 as well as the source position {9 s ^,(p s j). The observer’s 
collecting angle is chosen to be A 9 n = ±1°. To produce a statistically significant outcome 
even for observer’s polar angles 9 0 < 30° we consider a large number of sampling (a few 
lO 6 ) per source. However, the common feature they share is found to be a constant time 
separation of ~ 14 4/ between the major peaks of individual response functions, as found in 
FK08. To clearly illustrate these features of the response functions, we exhibit in Figure 5 
the response functions from flares at (a) 4> s ~ 4° and (b) ~ 180° for 9 S = 89° measured by an 
observer of 9„ — 60°. A constant time lag of ~ 144/ (equivalent to ~ 0.7 msec for a 104/. ) 
is clearly present. Due to multiple orbits around a fast-rotating black hole some photons 
make extra (integer-multiples) full orbits before reaching the observer. 

The response functions are different in the precise positions of the major peaks, as 
one would expect, but. as argued above, the time-lag between the major peaks is roughly 
constant. A t ~ 14 4/. regardless of the source location (9 s .o s ). The reason for this behavior 



was discussed in detail in FK08 (it remains the same in the present case too): it is due to 
the frame-dragging effects which force all photons emitted at any relative position between 
source and the observer, to reach the latter only by “swinging” around the black hole; in 
following such trajectories, a non-negligible fraction (~ 15 — 20%) of them reach the observer 
after multiple orbits around the black hole, inducing “light echoes” of nearly constant lag as 
shown in Figure 5. This is in agreement with similar relativistic calculations of broad iron 
line profiles by Beckwith & Done (2005) where they find that these higher-order photons 
contribute maximally ~ 20 — 60% of the total luminosity of the system. The persistent 
presence of such a constant time-lag in the response functions is the key factor leading to 
QPO features in the power spectra shown in §3. It should be noted that in the absence of 
strong frame-dragging (i.e. sources outside the ergosphere or Schwarzschild instead of Kerr 
black hole), the time interval among the peaks in the response functions is phase-dependent 
and precludes the presence of QPO in the power spectra for random source positions (see 
FK08). 

As discussed in §2.1 the X-ray interval between flares AT is Poisson-distributed with 
mean value comparable to the Keplerian period at the ISCO. Therefore, the collective re- 
sponse function from the ensemble of X-ray flares (i.e. the bolometric light curve) is obtained 
by summing over every contribution from all the sources at (8 Sti , 0 S|i ). That is, 

n 

= (18) 
1=1 

where n is a total number of random X-ray sources in our simulations. As an example, a 
synthetic light curve is shown in Figure 6a (i.e. superposition of all the response functions 
similar to Fig. 5) and its first 2 OAT msec segment in Figure 6b where AT = Af/(lOAf 0 ). 
In this specific case, the number of flares considered is n = 2000 with the source positions 
confined to near the equatorial plane (. 9 S = 89°) while the observer is set at 8 0 = 60°. Because 
of the random positions (<Tm) and bursting times of individual flares, the model light curve 
appears to be (in fact it is constructed to be) totally random with no apparent coherence 
(even in the absence of background noise). Nonetheless, the characteristic time-lag among 
the peaks in flux implicit in the form of the response function in Figure 5 is encoded in the 
light curve and it can be extracted (if present) by applying the standard time-series analysis 
techniques discussed in the next section §3. 


3. Power Spectra and Autocorrelation Functions 


In this section we discuss in detail the timing analysis metrics. ACFs and PSDs. of the 
model light curves produced using the prescriptions discussed above, pertaining to rapidly- 
rotating black holes: more specifically we study the dependence of the ACFs and PSDs on 
the sources' position and geometric structure as well as the position of the observer. 
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Time (M 1 sec) Time (IVh msec) 

Fig. 6. — (a) Model light curves constructed from the synthetic response functions similar 
to those in Figure 5. We assume n = 2000 random X-ray sources from 9 S = 89° for 9 a = 60°. 
(b) The first 20 M\ msec segment of (a). Values of the other parameters are the same as in 
Figure 4. 

As shown in §2 above, the form of the response functions and hence of the model light 
curves we produce depend mainly on the sources’ latitudes (given that we average over 
their azimuths Because the effects of frame-dragging are most prominent for sources 

near the equator, we first examine cases where most of the X-ray emission is concentrated 
preferentially near the equatorial plane, i.e. of 9 S)i = 89° (and corresponding positions in 
the lower hemisphere). Regarding the observer’s inclination angle 9 a we consider a wide 
range from nearly face-on ( 9 a ~ 30°) positions to nearly edge-on ( 9 a ~ 85°) ones. Since 
our primary goal in this investigation is to generalize the previous 2D results of FI<08, we 
restrict ourselves to a fast-rotating black hole case with a/M = 0.99 for which the source 
X-ray photons are subject to strong dragging of inertial frame (i.e. inside the ergosphere). 
Note that for a/M — 0.99 the ISCO radius is naturally well inside the ergosphere near the 
equatorial plane. 

We also consider, for a given observer’s angle 9 0 , the dependence of the ACF and PSD 
on the latitudinal distribution of the sources within the ergosphere. To this end we compute 
also the response function for source positions at different vertical heights with A 9 S = 9 ° by 
dividing the polar angle into five representative zones; i.e. 9 S = 89°, 80°, 71°, 62°, and 53° (see 
Fig. lb). We then examine the effects of the sources’ vertical extent by computing light curves 
from sources with incrementally larger vertical extent, i.e. from sources whose emission 
consists of the sum of point sources at a number of heights corresponding to the following set 
of inclination angles: 9 S — 9 S go (i.e. equatorial sources). (9 89 °4_ 80 o. 0s9 o +8o°+7i 0 - # 89 o + 80 °+ 7 i°+ 62 °. 
and fW+ 8 0 °+ 7 i°+ 62 o + 53 o - We do not consider sources at heights that would correspond to 
angles smaller than 53°, because for the given value of black hole spin ( a/M = 0 . 99 ) and 
a vertical source whose equatorial foot point is at the ISCO. a larger height would place 
that section of the source outside the ergosphere, thereby producing no contribution to the 
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QPO features via frame-dragging effect. Given that the effects of frame dragging decrease 
with increasing source height we expect the QPO effects to decrease with increasing vertical 
source extent. 

The timing properties of the model light curves are explored for different relative posi- 
tions between the source’s and observer’s polar angles by computing the ACF, R(t; 9 0 ,9 s ), 
and the corresponding PSD, P{u\ 0 O , 9 S ). For a discrete light curve I(t a , 9 a ; 9 S ) the normalized 
ACF is given by 


LJ 

R(ry,e„e,) = £ 


k~ 1 


I{tk,0o,0s)I(tk+j,0o,0 a ) 
P(t k ,9 o -0 S ) 


(19) 


where L is the number of time bins in the light curve of I. The corresponding PSD, in units 
of (rms/mean) 2 /Hz, is given by 


2T 

P(Vj\K9.) = 

7V ph 


(20) 


where T is the time duration of the light curve I and cij is the Fourier amplitude defined by 


1 L 

I(tk, 0 O ; 9,) = - ^ a 


j = i 


2nijk/L 


( 21 ) 


and iYph is the total photon counts in the entire light curve of /. Also, r, and Uj are timescale 
and frequency of a characteristic QPO signal, respectively. 


3.1. Dependence on Viewing Angle 

In this subsection we first present in Figure 7 a series of ACFs (left panels) and PSDs 
(right panels) as a function of the observer’s inclination (or viewing) angle 9 0 for model light 
curves from both near-equatorial sources, i.e. of 9 a — 89°, and a source that is vertically 
extended in height, modeled as a sum of individual sources at heights that correspond to 
polar angles 9 S = 89°. 80°, and 71°. These are shown in Figures 7a and b (equatorial source) 
and c and d (vertically extended source). 

As noted in FK08 the QPO signal in the light curves is imprinted by the frame dragging 
effects induced by the black hole spin. It is seen that a characteristic timescale of r ~ 
14 4/ = 0.734/i msec is persistently present for a wide range of inclination 0 O . which is a 
manifestation of the constant time-lag in the response function in Figure -5. Since these 
effects are more prominent for photon orbits that lie near the equator, this signature is more 
pronounced for arrangements that maximize the propagation of the received photons in this 
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Fig. 7. — Autocorrelation functions (left panels) and power spectral densities (right panels) 
for various observer’s inclination angles where 6 S = d 89 o (upper panels) and 6 a = # 89 o + 8 o o + 7 i o 
(lower panels). We show 6 0 = 30°, 40°. 50°, 60°, 70°, 80°, 85° from bottom to top curves with 
Au = 0.6 Hz bin size. Vertical positions of the curves are arbitrarily shifted for presentation 
purpose. Values of the other parameters are the same as in Figure 4. 


region. As shown in Figures 7a and b this effect is the most optimized for 0 O ~ 50° — 60°; 
for observers at smaller inclination (polar) angles the effect is weaker because the influence 
of frame-dragging is much reduced for photons that propagate close to the black hole spin 
axis. This becomes apparent by the decrease of the ACF peak at r ~ 0.73 Mi msec, or 
equivalently, the corresponding QPO amplitude with the decrease in 0 O (Fig. 7a and b). The 
QPO amplitude decreases also for 0 o > 60°; the reason for that is the increase of the fraction 
of the poloidal plane orbits that connect the source and the observer; these photons do not 
orbit around the black hole with coherence and as such they do not contribute to the QPO 
amplitude but in fact they dilute it. That is, the total signal is the result of a superposition 
of a mixture of both coherent and incoherent signals. 

Similarly, as shown in Figures 7c and d. a vertically extended source leads to reduction in 
the corresponding amplitudes of the ACF and QPO since the frame-dragging effects decrease 
as the source approaches the edge of the ergosphere. Comparison of Figures 7a and b with 
Figures 7c and d. which exhibit respectively the ACFs and PSDs for a set of observers at 
inclinations 0 o = 30°, 40°. 50°, 60°, 70°, 80°, 85° for an equatorial (a.b) and a vertically 
extended source (c.d); confirms these notions: the QPOs are more prominent for sources 
confined near the equator than for those which are vertically extended. The former also 
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persists over a wider range of values of 0 o , being strongest for 9 0 ~ 50° — 60°, while the 
latter, being intrinsically weaker, they disappear for 0 O > 60°. 

Finally, while the QPO and ACF secondary peak amplitudes do depend on the observer 
inclination, the corresponding QPO frequency and ACF peaks remain constant respectively 
equal to vqpo ~ 1. 4 M x 1 kFlz and T r ^ J 0.73A/i msec, a value representing the length of the 
photon paths around the ISCO, becoming considerably wider only for 0 O < 30°. 


3.2. Dependence on Source Concentration 

If the X-ray sources originate from processes other than those associated with the surface 
of a thin disk (e.g. flaring events due to magnetic reconnection, polar shocks or wind/jets), 
then the source may not be exclusively distributed in an equatorial fashion. In this subsection 
we examine in greater depth the effects of extending the source dimension vertically on the 
ACFs and QPOs respectively. To this end we have produced the light curves starting with 
an equatorially concentrated source (0 S = 89°) and consider the light curves from more 
extended sources by incorporating into the light curve also photons from sources at larger 
heights, or smaller values of the angle 9 S -, this is done incrementally by considering additional 
sources of the same intensity at the following set of source polar angles: (989°+80°, (989°+80 o +7i°i 
$ 89 o + 80 o + 7 i o + 62 °, and $ 89 ° + 80 ° +7i° + 62 ° +53° • We do not consider sources at any larger heights 
because for the conditions used here they would be outside the ergosphere and they would 
not contribute to the QPO signal by frame-dragging. 

In Figures 8a and b we present respectively the ACFs and PSDs for an observer at 
9 0 = 30° as the source height progressively increases to 9 S — 53° (i.e. sources extending all 
the way up to the static limit). Overall, the QPO signal does not appear to be very strong 
but it is present. We see that the additional (incoherent) photons dilute those with constant 
lag to cause the QPO to disappear with the source height. Similarly, in Figures 8c and d 
we present the ACFs and PSDs for an observer with high inclination of 9 0 = 60°. We see 
that in this case the QPO becomes again prominent when the X-ray sources are distributed 
near the equator as the photons that reach these directions are under a strong influence of 
frame- dragging. It appears that the latitudinal extension of the sources up to 6 S <71° can 
yield a discernible QPO signature, although it is relatively weak. 


4. Summary and Discussion 

In the present work we have extended the earlier 2D treatment of FK08 to consider the 
3D geometry of photon emission on the timing properties of accretion flows in the vicinity 
of rapidly spinning black holes. This has allowed us to explore the effects of placing both 
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Fig. 8. — Autocorrelation functions (left panels) and power spectral densities (right panels) 
for various source distributions where 9 0 = 6% o° (upper panels) and 9 0 = # 60 ° (lower panels). 
We show 9 S = 9ggo , $89 o +80°> #89 o +80 o +71° ! $ 89 o + 80 o + 71 o + 62 o > a nd $ 89 ° + 80 ° + 71 G +62° +53° from bot- 
tom to top curves with Au = 0.6 Hz bin size. Vertical positions of the curves are arbitrarily 
shifted for presentation purpose. Values of the other parameters are the same as in Figure 4. 


observers and sources at latitudes much different from the equatorial ones considered in FK08 
and, using these results, to consider also photon sources which are extended, i.e. sources 
whose vertical dimension is comparable or larger than their radial one. Our results are 
fundamentally consistent with those of FK08: we find that for the same reasons discussed in 
FK08, here too, frame-dragging causes the photons to follow trajectories that by and large 
lead to a significant fraction of them reaching the observer after an additional orbit around 
the black hole; as such, for a single X-ray flare, an observer sees multiple bunches of photons 
that arrive with a constant time-lag of order of ~ 144/ (i.e. light echoes ), independent 
of the source position relative to the observer. It is the independence of the lag on the 
source position that leads to the QPO in the signal of even a random ensemble of X-ray 
bursts, provided that they take place within the ergosphere. As elaborated in FK08. the 
resulting QPOs are different in character to those examined so far in the literature in that 
they do not require any modulation/oscillation as such in the signal. The difference in their 
character can be easily assessed in the data by calculating, in addition to the PSDs. also the 
ACFs, which in the present case has the form shown in Figures 7 and 8 . i.e. a peak near 
zero (self-coherence) with a second one at r ~ 144/ = 0.73x1/! msec (QPO), very different 
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from that of an oscillatory form, corresponding to QPOs due to some oscillation in the X- 
ray flux (Fukumura, Dong, Kazanas, & Shrader, in preparation). In summary, for sources 
concentrated near the equator (e.g. 0 S < 85°), we find the presence of persistent QPOs (and 
a series of harmonics) of uq PO ~ 1.4 kHz for a 10.1/., black hole (equivalently ~ 1.4 mHz for 
a 10'i\/@ AGN) for a wide range of line-of-sight viewing angles (from ~ 30° to ~ 80°). This 
QPO signal is the direct consequence of photon trajectories undergoing multiple rotations 
around the black hole due to its strong frame-dragging effects near the equator. As the 
source extends toward the mid-latitudinal region (~ 53°) the QPO signal is weakened by 
more dominant incoherent signals. 

The QPOs proposed in this work should be viewed as a new class of QPOs in view of their 
(i) expected frequency bandpass and (ii) their underlying mechanism, generic to the dragging 
of inertial frames on individual photons. As such, we do not believe these are associated 
with the observed HFQPOs often detected in accreting black hole systems (whose frequency 
is lower by a factor of ~ 5 and some of which exhibit a 2:3 frequency commensurability 
that the present model does not provide for) . The fundamental requirement for the presence 
of QPOs of the type discussed herein is that the X-ray flares would take place within the 
ergosphere of a rapidly-rotating black hole. The requirement that the ISCO lies within the 
ergosphere, then, implies that the black hole spin be a/M > 0.94. While this appears to be 
a rather strong constraint, one should point out that broad Fe line fits by emission from thin 
accretion disks around Kerr black holes, provide consistently estimates of the spin parameter 
near a/M ~ 0.99 (e.g. Brenneman & Reynolds 2006). This value is consistent with that used 
in the present work. 

The above results are relevant if the accreting gas surrounding the central engine is rela- 
tively tenuous and optically transparent to radiation from these X-ray sources. For example, 
this should be the case for optically-thin ADAFs where the accreting plasma distribution is 
quasi-spherical (though most of emission is still originating from the midplane of the gas), in 
contrast to the equatorial structure of a standard, thin accretion disk (e.g. Shakura & Sun- 
yaev 1973; Novikov & Thorne 1973). Such geometrically thick, optically thin configurations 
appear to be associated with objects such as our Galactic Center source Sgr A* (e.g. Yuan 
et al. 2003; Meyer et al. 2006), and Low Ionization Nuclear Emission-line Region sources or 
LINERs (e.g. NGC 4258; see Quataert 1999) and nearby giant ellipticals in low-luminosity 
AGNs or LLAGNs (e.g. M87; see Di Matteo et al. 2003). If the central emission region 
in these sources is indeed optically-thin characterized by ADAFs or its more generalized 
class of radiatively- inefficient accretion flows (RIAFs; see. e.g. Yuan et al. 2003 ). then one 
may expect a transparent environment in which photons can escape the production region 
without much additional scattering. 

QPO features have recently been associated with the PSD of AGN: X MM- Newton ob- 
servations of the the bright nearby Seyfert 1 galaxy Mrk 766. known to be a NLSl, exhibited 
statistically significant peak in the PSD of its light curve (Markowitz et al. 2007), modeled 


as Lorentzian of vqpo ~ 0.46 mHz, similar to QPOs detected from other black hole binary 
systems. Assuming an estimated mass of the Mrk 766 nucleus (~ 5 x lO'hY/©), suggested 
from optical H j3 emission line velocity dispersion (e.g. Wandel 2002), we can estimate the 
frame-dragging QPO frequency to be vqpo ~ 2.7 mHz, i.e. a factor of ~ 5 higher than the 
observed frequency. It is possible that other NLSls with slightly more massive nuclei (order 
of 1O 7 M 0 ) may exhibit a QPO frequency that can be explained in the context of our model. 
For example, a radio-loud NLS1, PKS 0558-504, shows a rapid X-ray flare presumably with 
a heavier nucleus of M ~ 4.5 x 10 7 A/ Q (e.g. Wang et al. 2001) which may make this source a 
good candidate for testing our prediction. Recently, Meyer et al. (2006) discussed a charac- 
teristic time-lag due to higher-order emission from our Galaxy Center (Sgr A*), which may 
also be a promising target for our investigation (because of larger mass and longer timescales 
compared to stellar-mass black holes). In addition to Seyfert Is above, perhaps intermediate 
mass black holes possibly associated with ULXs may also be promising sites to search for the 
QPOs proposed herein (e.g. Strohmayer et al. 2007). Although it may not be observationally 
easy to disentangle a potential QPO signal from high luminosity continuum (presumably at 
near-Eddington rate) from these objects, it remains to be studied. 

Our discussion so far has been confined to X-ray flare sources that orbit the black hole 
in the azimuthal direction with no significant radial motion. One could argue, however, that 
any source in the vicinity of an accreting black hole is likely to be plunging-in radially at a 
good fraction of c. This is certainly the case for sources interior to the ISCO. As pointed out 
above, in order that a source be both at the ISCO and within the black hole ergosphere, the 
spin parameter a should be a/M > 0.94. For more moderate values of the spin parameter, 
say, a/M ~ 0.9, the ISCO is outside the ergosphere and plasma traversing the latter should 
have considerable radial velocity in addition to its azimuthal motion. Since there is no 
reason for which this plasma could not produce X-ray flares (for example, shocked plasma in 
plunging regions can in principle serve as a good candidate for such a thermal/nonthermal 
X-ray sources, see Fukumura et al. 2007), the entire analysis described in this work can be 
extended to this circumstance too. It is expected, however, that because of its (substantial) 
radial velocity, the emitted photons would be Doppler beamed in the radial direction too and 
a large number of them would be lost through the black hole horizon; as a result, the QPO 
production efficiency may decrease to an unobservable level for sufficiently small values of 
the spin parameter a. A detailed study of this arrangement and the limiting value for which 
this approach can lead to QPO signals is deferred to future work. 

A related issue is that of the detectability of QPOs produced in the way described 
above. In this work we neglected the background distribution of signal/noise associated 
with other physical processes (e.g. emission from accreting flows and/or corona), photon 
counting statistics (Poisson noise) and also instrumental responses to the signal. Depending 
on the QPO signal strength relative to the externally contaminating (e.g. accreting gas 
and coronal) emission intensity, the induced light echo can be sufficiently weakened to a 


21 - 


statistically indiscernible level. To assess a potential observability of the predicted QPO 
signal in this scenario we need to combine the synthetic signal with the stochastic noise. If 
fraction of the observed X-rays does not participate in the production of the light echoes 
discussed herein, one should add to our synthetic signal an underlying constant (stable) 
component with some signal-to-noise ratio. Finally, in order to gauge the observability of 
the signal we propose, one should also add the noise associated with photon statistics and 
detector background. Such a detail treatment of the observability of the QPOs proposed in 
this work is deferred to another future work. 

We would like to thank the anonymous referee for a number of useful and insightful 
suggestions. K.F. is grateful to Fotis Gavriil for helpful discussions about timing analysis. 
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5. Appendix A 

In this Appendix we describe in details the notion of local isotropy and derive the 
equations necessary to calculate the subsequent photon trajectories in Kerr geometry. 

In order to correctly describe local isotropy in emission from a source orbiting at rela- 
tivistic speed, we first consider the differential photon distribution as a function of the angle 
ip between the source direction and the photon emission in the fluid frame 


dN d,N d ( cos ip) 

dil> d( cos ip) dip 


( 22 ) 


The notion of local isotropic emission is the requirement that the number of photons per 
local solid angle should be the same for all the direction of emission, i.e. that 


dN 

dQ 


No , 


(23) 


where No denotes the number of photons within the solid angle dfl which must be conserved 
from one frame to another and dfl = 2tt sin ip dip. Therefore, we obtain 


dN 

dip 


2ttAq sin ip , 


(24) 


which provides the weighting factor as a function of ip. 

In the LNRF (or lab frame), where the photon trajectories are computed, we similarly 


find 


dN dN f/(cos (.•) (/(cos 
dip' t/(cos u) d(cos ip') dip' 


(25) 


where again 


_dN_ 

(/(cos V 


= 2ttA 


a jV 0 


( 26 ) 


This preprint was prepared with tire AAS IATRX macros v5.2. 


( 27 ) 
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Using equation (9), we can express equation (25) as 


dN 

dip' 


2nN 0 f(p, ip) , 


where 


f(P, U) = 


(1 + (3 cos ip) sin ip 
(1 -/3 2 )V2 ’ 


( 28 ) 


Note that dN/dxp' — ► dN/dip as {3 — > 0 and we can see fj (dN j dtp') dip' = (dN / dip) dp’ = 

An N q regardless of the value of j3 as expected. Hence, f(fd,ip) describes the differential 
photon distribution per unit opening angle (i.e. the weighting factor) in general inertial 
frames. 


One can easily convert the above distribution into the differential photon distribution 
over a finite bin A ip' (in the LNRF) given by equation (10). Note that A ip' is not constant 
and one finds 


-y—^\ A ip' = A 'ib = 2 ttN 0 sin ip Nip 

d ip J \ dip J 


(29) 


which states that the number of photons emitted within the corresponding angular bin in 
the two frames is indeed the same (i.e. independent of inertial frame), as one expects. Then, 
it is guaranteed that for the entire sky 


"-E 




A ipi = 2ttN q I sin ipdip = 4 ttN 0 , 

Jo 


(30) 


where the integer i represents a discretised grid point in the polar opening angle and n 
denotes the total number of bins of the same angle from 0 to n. This prescription of photon 
distribution allows one to compute the number of photons in each local opening angle ip ' , 
consistent with isotropic emission at the fluid frame. For each such polar angle, we then 
launch photons in the azimuthal (about the source velocity) direction \ (see Fig. 2 ) between 
0 and 2 tt in equal intervals of width 


A\ = 2 ;r 



1 

Ao sin -ipAii> 


( 31 ) 


with this prescription, then, guaranteeing the proper number and photon directions in the 
LNRF consistent with isotropic emission in the fluid rest frame. 

Besides the Lorentz boost of the angles to the LNRF one must in addition relate these 
angles to the photon angular momenta (£, //), as these are the quantities needed in the 
integration of the photon geodesics. A simple geometry of this problem relates the opening 
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angle (,'/ and the azimuthal angle \ to the corresponding polar and azimuthal angles, S' and 
y. as described in §2.3 

cos •(// = sin S' sin 7' and cos S' = sin </7 cos \ , (32) 

which is schematically illustrated in Figure 2. 


6. Appendix B 

In Appendix A we discuss how to prescribe local isotropy from an arbitrary X-ray source 
position in Kerr geometry, which enables us to calculate the emission direction (5', 7') of an 
individual photon. In order to subsequently integrate the geodesic equations (photon trajec- 
tories) per photon given by equations (3) through (6) one needs to translate such a directional 
information into the photon’s angular momentum (equivalent to impact parameters) in the 
LNRF. Since we are dealing with 3D null geodesics each photon is characterized by two 
constants £ and r/ where is the axial component of angular momentum and 7/ is closely 
related to the polar component of angular momentum pa by 77 = pj + £ 2 cot 2 9 — a 2 cos 2 9 
(e.g. Bardeen et al. 1972; Chandrasekhar 1983). Here we derive analytic expressions in the 
most general case for (£, ?/) in terms of the local emission angles (S\ p r ). 

Considering equations (16) and (17) where the left-hand-side is related to the local 
emission angles and the right-hand-side to geometry and velocity, one can eliminate the 
derivatives ( i , r, 0. d>) using equations (3) through (6) to express the right-hand-side in terms 
of (£, rf). By solving these equations one can obtain analytic forms of £ = £(5', 7') and 
77 = p(S' , 7') explicitly as 

(1 — p 2 )ujA(sm 2 S' ± {(1 — p 2 )AS 2 C 2 sin 2 (5 / csc 2 7 / } 1/ “ 


Wr/) = 


V&8W) 


AS 2 esc 2 7' — (1 — p 2 )u; 2 A‘ 2 sin 2 S' 
p 2 CA{y - ci ~ ( 1 - p 2 )} + (1 - p 2 )(C - cm4£) 2 cos 2 <F 


(33) 


(34) 


(1-// 2 )CA 

where = i A + u 2 {p 2 A + r(r + 2M)} and p = cos 6 with all the quantities to be evaluated at 
the source position (r s ,0 s , <p H ). The plus/niinus sign for the second-term in the numerator of 
equation (33) depends on the initial direction of photons emitted in such a way that £ > 0 
for all emission; i.e., plus for photons initially emitted (locally) in the forward direction 
(0 < p' < 7T ) and minus for backward direction (~ < 7 ' < 2t r) with respect to the LXRF 
(which can be in rotation relative to the observer). Note that there is an "oifset" factor in 
directional information due to this rotation of local inertial frame [due to the quantity a.- in 
cqns. (33) and (34)]. In the absence of frame-dragging (a/M — * 0) we see that 

1/2 


c 

S 


±r 


1 




M/,p 


sin 6 sin - 


(35) 
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which is positive for both forward and backward emission with sin S' > 0, as expected (see 
Fig. 2). This also reproduces correct axial component £ in Newtonian geometry (r — > oc); 
when emission is exactly radial, either away (y 7 = 0) or toward ( 7 ' = 7r) the rotation 
(symmetry) axis, £ will vanish as required. When emission is completely parallel to the 
rotational axis (either S' = 0 or tt) , one also recovers £ — > 0. 

As a special case we can see that for emission from an 011 -axis (z-axis) source, i.e. as 
9 -► 0 (p -> 1) 

£ — > 0 and q — * finite , (36) 

indicating the absence of axial angular momentum, as anticipated. For another special case 
where emission is confined in the equatorial plane or (x,y)-plane as 9 — > 7t/2 (p — > 0) and 
S' — > 7t/2 we see 


£ — > finite and q —* 0 , 

where q, which in this case is exactly p#, vanishes, also as expected. 


(37) 



